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Abstract To investigate the effect of wall thickening 
around cell corners on the tangential Young’s modulus of 
coniferous early wood, tapered beam cell models in which 
the variation of the cell wall thickness in the axial direction 
was taken into account were constructed for seven species. 
Their tangential Young’s moduli were compared with the 
experimental results. The calculated Young’s moduli of 
tapered beam cell models were larger than those of the 
models composed of the cell walls with uniform thickness, 
although both models showed almost the same density. For 
some species the calculated Young’s moduli of the models 
in which the cell wall thickness increased curvilinearly in the 
axial direction were much closer to the experimental values. 
The reduction of the radial cell wall deflection due to the 
increase of the stiffness around cell corners was considered 
to increase the tangential Young’s modulus of a wood cell. 

Key words Tangential Young’s modulus ■ Wood cell model 
■ Cell wall thickness • Cell corner 


Introduction 

In previous articles 1 ' 2 the relation between the tangential 
Young’s moduli and the transverse cell shapes for seven 
kinds of coniferous early wood was investigated using cell 
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models constructed by power spectrum analysis. The calcu¬ 
lated Young’s moduli of the cell models explained qualita¬ 
tively the change of the experimental Young’s moduli 
with density and the difference in the experimental values 
among species. However, the calculated values were lower 
than the experimental values for all species exa m ined. Two 
reasons may be considered for the lower calculated values 
of the tangential Young’s modulus of the cell model: the 
validity of the axial and bending Young’s moduli of the 
wood cell wall in the perimetric direction and the variation 
of the cell wall thickness in the perimetric direction. 

The axial and bending Young’s moduli were calculated 
as 14.3 and 15.6GPa, respectively, in previous articles. 1,2 If 
these two values are about 20.0 GPa, the calculated Young’s 
moduli become larger and closer to the experimental re¬ 
sults. However, this value cannot be obtained from a cell 
wall model constructed with reference to a standard wood 
cell wall structure. 1 " 6 Therefore, the two values for Young’s 
moduli of the cell wall cannot be the main reason for the 
lower calculated values of the tangential Young’s modulus 
of the cell model. 

The cell models described in the previous articles were 
constructed using cell walls with uniform thickness in the 
perimetric direction. In reality, the cell wall thickness is not 
uniform, being somewhat thicker around cell comers. Tak¬ 
ing into account the effects of beam thickness on stiffness, it 
is estimated that the thickening of the cell walls around cell 
corners contributes to an increase in transverse Young’s 
modulus of a wood cell. The variation of the cell wall thick¬ 
ness in the perimetric direction may become an important 
factor when analyzing the mechanical behavior of a wood 
cell. 

In this study we derived an equation to calculate the 
tangential Young’s modulus of a tapered beam cell model in 
which the cell wall thickness varied in the axial direction. 
Furthermore, two types of tapered beam cell model were 
constructed for the species used in the previous articles, and 
their tangential Young’s moduli were compared with the 
experimental values to investigate the effects of tapered cell 
walls on the tangential Young’s modulus of coniferous early 
woods. 
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Fig. 1. Transverse shape of A-model. R, axial length of the radial cell 
wall; T, axial length of the tangential cell wall; 6 , element angle; f„, half 
of the cell wall thickness at the central portion; f c , half of the cell wall 
thickness at the cell corner 


Tangential Young’s modulus of a tapered beam 
cell model 

Figure 1 shows a tapered beam cell model and its geometric 
parameters. R and T show the axial length of radial cell wall 
and that of tangential cell wall, f 0 and t c show the cell wall 
thickness in the central portion and that at cell corner, and 
9 shows the element angle, respectively. When this model is 
loaded in the tangential direction, the radial cell walls bend, 
and all of the cell walls deform in the axial direction. 7 The 
deflection of the radial cell wall <5 is given by 


where b is the cell wall depth, and t R (x) is a function express¬ 
ing the variation of the radial cell wall thickness along the 
axial direction. The expression 2f R (x) is equal to the radial 
cell wall thickness at x. The tangential strain of the cell 
model caused by the radial cell wall deflection e b can be 
derived as 
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where o T is the stress in the tangential direction. 

The axial displacements of the radial and tangential cell 
walls, elfRsin# and ejT, respectively, are 
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where e a R and ej are the axial strain of the radial and the 
tangential cell walls, respectively, E a is the axial Young’s 
modulus of the cell wall, and t T (x) is a function expressing 
the variation of the tangential cell wall thickness in the axial 
direction. The expression 2 t T (x) is equal to the tangential 
cell wall thickness at x. The tangential strain of the cell 
model caused by the axial deformation of the cell walls e a 
can be derived as 
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All strain of tapered beam cell model is the sum of Eqs. (3) 
and (5). Therefore, the tangential Young’s modulus of the 
cell model E T is given by 
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where W, W , £ b , and I{x) are a force acting in the tangential 
direction, the virtual force, the bending Young’s modulus of 
the cell wall, and the second moment of inertia of the cell 
wall, respectively. Substituting W = 1 for Eq. (1) gives 
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(2) As a simple case, an equation for E T of a tapered beam cell 
model in which the cell wall thickness increased linearly 
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Table 1. Geometric parameters of A-models for seven species 


Species 

R {um) 

T (pm) 

0(°) 

t 0 (ftm) 

4 Om) 

t (pm) 

Cryptomeria japonica 

26.4 

19.1 

17 

0.99 

1.21 

2.20 

Chamaecyparis obtusa 

20.0 

18.7 

21 

2.03 

2.48 

4.50 

Picea glehnii 

20.2 

23.0 

13 

1.58 

1.93 

3.50 

Pinus densiflora 

25.1 

31.0 

12 

2.39 

2.92 

5.30 

Pinus radiata 

25.6 

28.0 

18 

2.48 

3.03 

5.50 

Metasequoia glyptostroboides 

28.7 

33.1 

10 

1.26 

1.54 

2.80 

Tsuga heterophylla 

22.2 

34.2 

6 

1.80 

2.20 

4.00 


R, axial length of the radial cell wall; T, axial length of the tangential cell wall; 9 , element angle; 
to, half of the cell wall thickness at the central portion; t a half of the cell wall thickness at the cell 
corner; t , cell wall thickness of S-model 1,2 


from 2f 0 to 2 t c (f 0 < f c ) in the axial direction, as shown in 
Fig. 1, was obtained as follows; 
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Hereafter, the cell model in Fig. 1 is referred to as A-model. 
£ t s of A-models for seven species (Chamaecyparis 
obtusa Endl., Cryptomeria japonica D. Don, Metasequoia 
glyptostroboides Hu et Cheng, Picea glehnii Mast., Pinus 
densiflora Sieb. et Zucc., Pinus radiata D. Don, Tsuga 
heterophylla Sarg.) were calculated using the geometric pa¬ 
rameters shown in Table l. 1,2 The parameters f 0 and t c were 
determined by two requirements: the relation of 

kZA = o.l (8) 

t 



where t is the cell wall thickness of the models constructed 
in the previous articles 1,2 (hereafter that model is referred to 
as S-model), and that the A-models showed the same den¬ 
sity values as those of the S-models. The deflection of lami¬ 
nates such as the wood cell walls significantly depends 
on the volume fractions and mechanical properties of 
the outer layers. Therefore, the elastic modulus of the cell 
wall calculated using the lamination theory 8 is considered 
to be proper for E b . Regarding the mechanical property 
of the cell wall as uniform along the axial direction, 
the values of 14.3 and 15.6GPa were used for E a and E b , 
respectively. 1,2 


Tangential Young’s modulus of B-model 

The thickness of real wood cell wall increases curvilinearly 
from the central portion of the cell wall to the cell corner, as 
shown in Fig. 2. An equation for E y of a tapered beam cell 
model in which this cell wall shape was taken into account. 


Fig. 2. Transverse cell shape of Pinus densiflora 


as shown in Fig. 3, was derived. The following equations 
were adopted for t R (x ) and f x (x). 
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Hereafter, the cell model in Fig. 3 is called B-model. Sub¬ 
stituting Eq. (9) for Eq. (6), the E r of B-model is given by 
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Fig. 3. Transverse shape of B-model 


Table 2. Coefficients of f R (.x) and c T (x) for seven species 


Species 

a 

Pr 

Pi 

m 

Cryptomeria japonica 

0.20 

18.0 

13.0 

0.50 

Chamaecyparis obtusa 

1.20 

12.2 

11.4 

0.50 

Picea glehnii 

0.85 

13.9 

15.7 

0.50 

Pinus densiflora 

1.60 

15.3 

18.9 

0.50 

Pinus radiata 

1.65 

15.2 

16.6 

0.50 

Metasequoia glyptostroboides 

0.50 

19.6 

22.6 

0.50 

Tsuga heterophylla 

1.00 

13.9 

21.5 

0.50 


Table 2 shows coefficient values of a, f} R , /? x , and m for 
calculating the £ x s of B-models for seven species. The den¬ 
sities of B-models calculated from these coefficients were 
almost the same as those of S-models. The values of 14.3 
and 15.6GPa were used for E. t and E b . 12 


Finite-element analysis of deflection and axial 
deformation of tapered cell walls 

The deflection and axial deformation of tapered cell walls 
used in A-models and B-models for seven species were 
investigated using a p-version finite-element analysis pro¬ 
gram (StressCheck V4.0; Engineering Software Research & 
Development.) and compared with those obtained from 
Eqs. (2) and (4). Figure 4 shows an element division of the 
tapered cell wall used in this analysis. X and Y indicate the 
axial and thickness directions of the cell wall, respectively. 



Fig. 4. Element division of tapered cel! wall for p-version finite- 
element analysis. X , axial direction; Y, thickness direction; P, point 
force; A, axial force 


All of the tapered cell walls examined were divided into five 
equal parts in the X direction. The cell wall was completely 
constrained at the end with the largest thickness (the right 
end in Fig. 4), and a point force in Y direction or an axial 
force in X direction was applied at the opposite end. A 
value of 14.3 GPa was used for the Young’s modulus E in 
this analysis. For Poisson’s ratio, a value of 0.3 was used. 
Furthermore, the £ x s of A-models and B-models based on 
the results of finite-element analyses were calculated and 
compared with those obtained by Eqs. (7) and (10). 


Comparison of experimental and calculated values of 
tangential Young’s modulus 

Figure 5 shows the relation between densities and experi¬ 
mental Young’s moduli and the calculated moduli of A- 
models and S-models for seven species. The difference was 
within 3.8% between cell wall deflections of A-models cal¬ 
culated from Eq. (2) and those obtained from finite-element 
analyses. For axial displacement, the difference was within 
1.7% between the values calculated from Eq. (4) and those 
obtained from finite-element analyses. The E T s of A-models 
for seven species are represented by the values calculated 
by Eq. (7), as they are almost the same as the results ob¬ 
tained by finite-element analyses. For all species examined, 
the E r of the A-model was larger than that of the S-model, 
although both models showed the same density. The £ T of 
S-model can be given by the following equation. 1 ' 27 

1 cosO(2T/R + sin : 0j j cos’ 6 

E t {t/R){T/R + sin G)E a ( t/R.f(T/R + sin d)E b 

( 11 ) 

The first terms of Eqs. (7) and (11), A’ and A, express the 
axial contribution; and the second terms, B 1 and B, express 
the bending contribution, respectively. A '/A and B'/B were 
calculated for seven species to investigate the effect of the 
tapered cell walls on the deformation of cells. Results are 











113 


Fig. 5. Comparison of experimen¬ 
tal and calculated values of tan¬ 
gential Young’s modulus. Squares, 
Cryptomeria japonica', circles, 
Chamaecyparis obtusa', diamonds, 
Picea glehnii', open crosses, Pinus 
densiflora ; triangles, Pinus radiate, 
inverted triangles. Metasequoia 
glyptostroboides ; squares with x, 
Tsuga heterophylla; plus sign, cal¬ 
culated value of S-model; x, cal¬ 
culated value of A-model 
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Table 3. Calculated values of A'/A and B'/B for seven species 


Species 

A'/A 

B'/B 

Cryptomeria japonica 

1.00 

0.87 

Chamaecyparis obtusa 

1.00 

0.87 

Picea glehnii 

1.00 

0.87 

Pinus densiflora 

1.00 

0.87 

Pinus radiata 

1.00 

0.87 

Metasequoia glyptostroboides 

1.00 

0.87 

Tsuga heterophylla 

1.00 

0.87 


A'/A, ratio of the first term of Eq. (7) to that of Eq. (11); B'/B, ratio of 
the second term of Eq. (7) to that of Eq. (11) 


shown in Table 3. Because A'/A shows a value of about 1.0 
for all species, the axial contribution is not influenced by 
the change from the straight cell walls to the tapered ones. 
On the other hand, B’/B shows a value below 1.0 for all 
species. These results indicate that a smaller deflection of 
the tapered cell wall increases £ x of cell models. 

Figure 6 shows the relation between densities and ex¬ 
perimental Young’s moduli as well as calculated moduli of 
B-models, A-models, and S-models for seven species. The 
difference was 3.4%-10.0% between cell wall deflections of 
B-models calculated from Eq. (2) and those obtained from 
finite-element analyses. For axial displacement, the differ¬ 
ence was within 1.9% between the values calculated from 
Eq. (4) and those obtained from finite-element analyses. 
Therefore, the E T of the B-model calculated from Eq. (10) 
and that obtained from finite-element analyses are indi¬ 
cated in Fig. 6 for each species. For all species examined, the 
E t of the B-model was the largest among the three cell 
models although densities were of almost the same value. 
The values of Chamaecyparis, obtusa, Pinus radiata, Picea 
glehnii, and Tsuga heterophylla became considerably closer 
to their experimental values. The t 0 and t c values of B- 
models calculated by Eq. (9) are shown in Table 4. The t 0 


Table 4. Half of the cell wall thickness of B-models for seven species 
calculated from Eq. (9) 


Species 

t 0 (wm) 

4 C«m) 

Radial wall 

Tangential 

wall 

Cryptomeria japonica 

0.70 

2.16 

2.11 

Chamaecyparis obtusa 

1.70 

3.97 

4.05 

Picea glehnii 

1.35 

2.68 

2.71 

Pinus densiflora 

2.10 

4.43 

4.38 

Pinus radiata 

2.15 

4.82 

4.84 

Metasequoia glyptostroboides 

1.00 

2.38 

2.35 

Tsuga heterophylla 

1.50 

3.48 

3.44 


t 0 , half of the cell wall thickness at the central portion; t c , half of the cell 
wall thickness at the cell corner 


Table 5. Calculated results of A'/A and B"/B for seven species 


Species 

A'/A 

£"/B Cal 

B"/B n 

Cryptomeria japonica 

1.09 

0.69 

0.68 

Chamaecyparis obtusa 

1.05 

0.74 

0.79 

Picea glehnii 

1.04 

0.74 

0.79 

Pinus densiflora 

1.04 

0.75 

0.83 

Pinus radiata 

1.04 

0.75 

0.82 

Metasequoia glyptostroboides 

1.06 

0.70 

0.69 

Tsuga heterophylla 

1.05 

0.72 

0.74 


A'/A, ratio of the first term of Eq. (10) to that of Eq. (11); B"/B, ratio 
of the second term of Eq. (10) to that of Eq. (11); Cal, calculated from 
Eqs. (10) and (11); FE, calculated by finite-element analysis 


value of the B-model was smaller than that of the A-model 
and the f/2 for each species. Table 5 shows the calculated 
results of A"/A and B"/B for each species where A" and B" 
express the first and second terms of Eq. (10), respectively. 
The value of A"/A was slightly larger than that of A'/A for 
all species. Flowever, this increase may not significantly 
influence the deformation of a cell because A is a smaller 
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Fig. 6. Comparison of experimen¬ 
tal and calculated values of tan¬ 
gential Young’s modulus. Stars, 
calculated value of B-model; Cal, 
calculated from Eq. (10); FE, cal¬ 
culated by finite-element analysis. 
See Fig. 5 for explanation of other 
symbols 


cS 

Ph 

o 

05 

3 

3 

-o 

o 


05 

'00 

s 

3 

O 

3 

o> 

00 

3 

03 

H 


0. 

0 . 

0. 

0 . 

0. 



Q I _ _,_I_l 

0.0 0 . 1 0.2 0.3 0.4 


Density y (g/cm 3 ) 


value than B. 1 ' 2 On the other hand, the value of B"/B was 
smaller than that of B'/B for all species; that is, the cell wall 
deflection of the B-model decreased compared with that of 
the A-model although t 0 of the B-model was a smaller value. 
From these results, it is considered that the larger values of 
£ t of the B-models are mainly attributed to the decrease in 
cell wall deflection. 

When a wood cell is loaded in the tangential direction, 
the radial cell walls deflect owing to the displacement of cell 
corners. The wall thickening around cell corners increases 
the stiffness and reduces the displacement, resulting in a 
reduction of bending deflection of whole radial walls. At the 
same time, the cell walls deform in the axial direction. FIow- 
ever, the axial contribution to cell deformation is negligible 
because the whole strain of a cell by this deformation is 
much smaller and the wall thickening around cell corners 
does not influence significantly the axial deformation of the 
cell wall. Therefore, in coniferous early woods, the effect of 
the large stiffness around cell corners on the bending prop¬ 
erty of the radial cell wall is considered to relate strongly to 
the larger experimental value of tangential Young’s modu¬ 
lus than the calculated value of the S-model. These consid¬ 
erations suggest that the wall thickness around cell comers 


as well as transverse cell shapes are important factors in the 
mechanical properties on the transverse plane in wood. 
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